COMPARATIVE ANALYSIS FOR PAIR OF DYNAMICAL 
SYSTEMS, ONE OF WHICH IS LAGRANGIAN. 



R. A. Sharipov 

Abstract. It is known that some equations of differential geometry are derived 
from variational principle in form of Euler-Lagrange equations. The equations of 
geodesic flow in Riemannian geometry is an example. Conversely, having Lagrangian 
dynamical system in a manifold, one can consider it as geometric equipment of this 
manifold. Then properties of other dynamical systems can be studied relatively as 
compared to this Lagrangian one. This gives fruitful analogies for generalization. In 
present paper theory of normal shift of hypersurfaces is generalized from Riemannian 
geometry to the geometry determined by Lagrangian dynamical system. Both weak 
and additional normality equations for this case are derived. 



1. Introduction. 

Let M be some smooth manifold and let L be Lagrange function of some La- 
grangian dynamical system in M. This is smooth scalar function depending on 
point p of M and on vector v at this point. In other words, L is a function of point 
q = (p, v) of tangent bundle TM. We treat L as basic equipment of manifold M 
(like metric tensor in Riemannian geometry or simplectic structure in simplectic 
geometry). Therefore we assume L to be satisfying some special requirements. Let 
x , . . . , x n be coordinates of point p in some local chart in M, and let v , ... , v n 
be components of vector v in this local chart: 

Lagrange function L in local coordinates x 1 , ... , x n is represented as the function 
of 2n arguments L — L{x , , . . . , x n ,v , . . . , v n ), while corresponding Lagrangian 
dynamical system is given by the following ODE's: 

dt\dv l ) dx l y ' 

Due to x % = v l in (1.2) vector (1.1) is called velocity vector. In addition to (1.2) we 
consider so called modified Lagrange equations: 



v 



d fdL\ _ 1 dL 
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Denominator fi in these equations (1.3) is determined by the following formula: 

« = !>•§■ m 

i— 1 

Like original Lagrange equations (1.2), modified equations (1.3) with denominator 
(1.4) arise in applications. As shown in [1], they describe wave front dynamics for 
various wave propagation phenomena. Now let's denote by p the covector at the 
point p € M given by its components in local chart: 

Covector p is called momentum covector. Pair composed by point p and covector 
p is a point of cotangent bundle T*M. Therefore (1.5) defines a map 

X:TM^T*M. (1.6) 

This map is called Legendre transformation. It is well known in mechanics (see [2]). 

Note that first equations in (1.2) and (1.3) are different. Therefore vector v 
cannot be interpreted as velocity vector for modified dynamical system (1.3). For 
this purpose we introduce vector u with components 

v i 

W = Ty (1.7) 

Formula (1.7) determines another map, which is similar to (1.6): 

H : TM — > TM. (1.8) 

Definition 1.1. Lagrangian dynamical systems (1.2) and (1.3) and their Lagrange 
function L are called regular if both maps (1.4) and (1.8) are diffeomorphisms and 
if denominator f2 in modified Lagrange equations (1.3) given by formula (1.4) is 
positive at all points q — (p, v) of tangent bundle TM, where |v| ^ 0. 

In order to get simplest example of regular Lagrangian dynamical system one 
should assume M to be equipped with Riemannian metric g and one should take 
Lagrange function L to be quadratic with respect to velocity vector v: 

-. n n 

i^EE^v-tv,...,*"). (i.9) 
»=i j=i 

Though it is simple, this example covers all real mechanical systems if we neglect 
friction in junctions and other forms of energy pumping and dissipation. First term 
in (1.9) is kinetic energy, while second term U(x 1 , . . . ,x n ) is potential energy of 
mechanical system. For L of the form (1.9) Legendre transformation (1.6) looks 
like index lowering procedure in metric g. Indeed, applying (1.5) to (1.9), we find 

n 

I'. ( L1 °) 
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If L is treated as geometric equipment of manifold M, then for L of the form (1.9) 
this equipment is equivalent to Ricmannian metric in essential. This is the case 
considered in papers [3] and [4]. In present paper we consider more general case, 
when geometric equipment of manifold M is determined by some arbitrary regular 
Lagrange function L, which is not necessarily given by formula (1.9). 

Now let's consider another dynamical system. It can be either Lagrangian or not 
Lagrangian, but we assume it to be second order dynamical system. More precisely, 
we assume that it is given by differential equations 

x l =u\ ii l = ^(x 1 ,... ,x n ,u 1 ,... ,u n ). (1.11) 

We call (1-11) Newtonian dynamical system, since these equation are similar to 
those expressing Newton's second law for the motion of a particle of unit mass 
m = 1 under the action of force Certainly, we cannot touch all problems 
associated with pairs of dynamical systems of the form (1.2) and (1.11). In present 
paper we construct theory of normal shift for such pairs of dynamical systems, thus 
realizing the project claimed in previous paper [1]. 

Theory of Newtonian dynamical systems admitting normal shift was initiated 
in 1993 in preprint [5]. At first dynamical systems (1.11) in Euclidean space M™ 
were considered, then theory was extended for dynamical systems in Riemannian 
and Finslerian manifolds. This phase, which lasted 7 years from 1993 till 1998, is 
reflected in theses [6] and [7] (see also appropriate references therein). For recent 
results in the theory of dynamical systems admitting normal shift see papers [8-12] 
and papers [1], [3], and [4] already mentioned above. Below in section 2 we give 
some preliminary information and motivated definitions. Then in further sections of 
this paper we construct theory of Newtonian dynamical systems (1.11) admitting 
normal shift for the case of manifolds, geometry of which is not Euclidean, not 
Riemannian, and even not Finslerian, but is given by Lagrange function L of some 
regular Lagrangian dynamical system in them. 

2. Normal shift of hypersurfaces. 

Suppose that S is some arbitrary smooth hypersurface in M. We say that S is 
equipped with smooth transversal vector field if at each point p G S some nonzero 
vector u(p) transversal to S is fixed. In this case one can consider the following 
initial data for Newtonian dynamical system (1.11): 

ac*| =x l ( P ), uH =u l (p). (2.1) 

lt=o lt=o 

Here x l (p) are coordinates of point p in some local chart in M and u l (p) are compo- 
nents of transversal vector u(p) in this chart. Applying initial data (2.1) to (1.11), 
we get a family of trajectories of this dynamical system starting at the points of S. 
In local chart it is represented by functions 

' x 1 = x 1 ^^), 

< (2-2) 

^ x n = x n (t,p). 
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These trajectories are transversal to S. If we fix time instant t ^ and gather all 
points of trajectories (2.2) corresponding to this time instant (see Fig. 2.1), we get 

another hypersurface S t and diffcomorphism 

ff-S^St (2.3) 

binding St with initial hypersurface 1 . Diffco- 
morphism ft (or, more precisely, the whole 
set of diffeomorphisms f t ) is called a shift of 
S along trajectories of Newtonian dynamical 
system (1.11). Note that the shift (2.3) keeps 
transversality in local. This means that trajec- 
tories of shift are transversal not only to initial 
hypersurface S, but to all hypersurfaces S t for 
sufficiently small values of t. Shift f t is called 
normal shift if it keeps orthogonality of St and 
trajectories in some sense. In previous papers 
(see [8-12] and earlier) orthogonality was understood in the sense of some metric 
either Euclidean, Ricmannian, or Finslcrian. Below we shall understand it in the 
sense of Lagrange function L as it was suggested in paper [1]. 

Let y 1 , ... , y™ -1 be some local coordinates on initial hypersurface S. Due to 
diffeomorphisms of shift (2.3) they can be transferred to all hypersurfaces St- Func- 
tions (2.2) in terms of local coordinates y 1 , . . . , y"" 1 are written as follows: 



f x 1 =x 1 (t,y\... ,y n - 1 ), 



{ x n = x n (t,y\... ,y n - 1 ). 



(2.4) 



Their time derivatives u % — x % are components of velocity vector u. It's easy 
to understand that partial derivatives of these functions with respect to y % are 
components of vector tangent to St- It is called z-th vector of variation (or, more 
precisely, i-th vector of variation of trajectories). We denote this vector by Tf. 

^=T}g^ + --- + r?g^, where T t = — (2.5) 

Vectors n, ... , t„_i form a base in tangent space to S t - For normal shift they 
should be perpendicular to velocity vector u. According to receipt from paper [1], 
we define the following deviation functions 



= ( P \T t )=J2PsT l S . (2.6) 



Here we have no metric, therefore scalar product (p | r^) is nothing else, but sym- 
bolic notation for the sum in right hand side of (2.6). 



1 In general, this is true only for sufficiently small hypersurface S and for time instants t 
sufficiently close to initial time instant t = 0. However, our further consideration is local. Therefore 
here we shall not discuss the problem of globalization for diffeomorphisms (2.3) referring reader 
to paper [12], where some aspects of this problem are studied 
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Definition 2.1. Shift of hypersurface S along trajectories of Newtonian dynamical 
system (1.11) denned by initial data (2.1) is called normal shift if all deviation 
functions ifi given by formula (2.6) are identically zero. 

Deviation functions <pi, ... , ip n -i are used as a measure of deviation of shift ft 
from being a normal shift. Their vanishing is indicator of normality. 

3. Relative form of the equations of Newtonian dynamics. 

Note that covector p in (2.6), according to (1.5), depend on components of vector 
v. However, v is not velocity vector for dynamical system (1.11). Vectors v and 
u are bound by the map (1.8), which is expressed by formula (1.7) in local chart. 
This map is diffcomorphism due to our assumptions (see definition 1.1). Therefore 
we can transform (1.11) to variables x 1 , ... , x n , v 1 , ... , v n . Here we write 



v l = tf^jc 1 ,.-- ,x n ,v\... ,v n ). (3.1) 



One could express "J/ 1 , ... , \E ,n through functions $ 1 , ... , $" in (1.11). However, 
the latter ones are arbitrary functions. Therefore we can assume vp 1 , . . . , in 

(3.1) to be arbitrary functions as well, with no need to follow their relations to the 
functions <i> 1 , . . . , $™ in (1.11). 

In the next step we use Legendre transformation (1.6) in order to transform 
differential equations (3.1) further. Now we write them as 

±l = n' Jt (cv ) " n = Qtl (3 - 2) 

where Qi — Qi(x 1 , . . . ,x n 1 v 1 , . . . ,v n ). Since A is diffcomorphism, these equation 

(3.2) are equivalent to previous ones (3.1). This is relative form of the equations 
of Newtonian dynamics (1.11). Relative, because in writing them we need another 
dynamical system (1.3). 

4. Extended tensor fields. 

Let's consider quantities Qi, ... , Q n in (3.2). By means of direct calculations 
one can check up that these quantities are transformed as components of covector 
under the changes of local charts in M. They define a covector at the point p, 
where p is a point with local coordinates x , ... , x n . Let's denote this covector by 
Q and note that it depends not only on p, but on components of vector v as well. 
This means that Q fits the definition of extended covector field. 

Definition 4.1. Extended covector field X on a manifold M is a covector- valued 
function that to each point q = (p, v) of tangent bundle TM puts into correspon- 
dence some covector from cotangent space T*{M) at the point p = Tr(q) in M. 

Similarly one can define extended tensor fields on the manifold M. First let's 
consider the following tensor product of tangent and cotangent spaces: 



7J(p, M) = T p {M) ® . . . ® T p (M) ® T;{M) ® . . . ® T p *(M) . 
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Tensor product TJ(p, M) is known as a space of (r, s)-tensors at the point p of the 
manifold M. Pair of integer numbers (r, s) determines the type of tensors. Ele- 
ments of the space TJ(p, M) are called r-times contravariant and s-times covariant 
tensors, or tensors of the type (r, s), or, for brevity (r, s)-tensors. 

Definition 4.2. Extended tensor field X of the type (r, s) on a manifold M is a 
tensor- valued function that to each point q = (p, v) of tangent bundle TM puts 
into correspondence some tensor from tensor space TJ(p, M). 

As far as I know, extended tensor fields first arose in Finslerian geometry. They 
was intensively used in theory of Newtonian dynamical systems admitting normal 
shift (see theses [6], [7] and references therein). Their application to Lagrangian 
and Hamiltonian dynamical systems is explained in [1] , [3] , and in [4] . 

5. Momentum representation for extended tensor fields. 

Note that we can replace tangent bundle TM in definition 4.2 by cotangent 
bundle T*M. Then we obtain another definition of extended tensor field. 

Definition 5.1. Extended tensor field Y of the type (r, s) on a manifold M is a 
tensor- valued function that to each point q = (p, p) of cotangent bundle T*M puts 
into correspondence some tensor from tensor space TJ(p, M). 

Extended tensor fields as given by definitions 4.2 and 5.1 are two different objects. 
However, they can be related to each other by Legendre map (1.6) that links tangent 
bundle TM and cotangent bundle T*M. Suppose that 

X = YoA, Y = X°A" 1 . (5.1) 

If relationships (5.1) hold, we say that Y is p-representation (or momentum 
representation) of extended tensor field X. Similarly, in this case we say that X is 
v-representation (or velocity representation) of extended tensor field Y. 

6. Weak normality condition. 

Let's consider Newtonian dynamical system written in relative form (3.2) and 
let's consider some one-parametric family of trajectories p = p(t, y) of this dynam- 
ical system. In local chart this family of trajectories is represented by functions 

x 1 = x^y), 

(6-1) 

x n = x n (t,y) 

(compare with (2.4)). Time derivatives of (6.1) define vector v according to first 
part of the equations (3.2). Components of v depend on t and y: 

v 1 =v\t,y), v n = v n (t,y). (6.2) 

Like in (2.5), we can define variation vector r. It is given by formula 



(6.3) 
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Here we have only one variation vector since, besides time variable t, in (6.1) we 
have only one parameter y. Therefore we have only one deviation function 

n 

^=(p\r)=Y / Psr s . (6.4) 

8=1 

Let's differentiate (6.2) with respect to parameter y. This yields a series of functions 

9 i = 9 i (t,y) = ^. (6.5) 

Double set of functions r 1 , . . . , r" and 6 1 , ... , 9 n satisfy a system of linear ordi- 
nary differential equations with respect to time variable t. This system is obtained 
as linearization for (3.2). Differentiating first part of the equations (3.2) with re- 
spect to parameter y, we get the following expression for time derivative f % : 

n ^n 2 dv s 

s=l 

(6.6) 

y, y, d 2 L V 1 V k 6 S yy y> d 2 L V 1 V k T 8 

fc=l s=l fe=l s=l 

Differentiating second part of the equations (3.2) with respect to y, we obtain 
™ f) 2 r " F) 2 t ™ " f) 3 r 

y o s + y f s + yy v k e s + 

^— J dv l dv s t-* 1 dv l dx s ^— J dv l dv s dv k 

s=l s=l k=l s=l 

+ yy Jii: ± k 9 s + y y . « fc r s + 

Z_^l Q y l Q y s Q x k Z_^i ^2 Qyk 

k=l s=l k=l s=l 

y,y, d 3 L . k _ s y, d 2 L 6* y, Q 2 L T° 

+ Z^l^Q v i Q x s Q x k X T 1^ Q x i Q v s n 1^ Q x i Q x s ^ + (6.7) 
k—1 s—1 s—1 s—1 

y, dL dL 9 s y, y, dL d 2 L v k 9 s 

s=l k=l s=l 

Tyy, PL 8 2 L V k T S = y^8Qi a 9Qi gs 

^ ^ dx l dv k dx s Q 2 ^ dx s T ^ dv s 

k—1 s — 1 s—1 s—1 

The equalities (6.7) are not resolved with respect to derivatives 9 1 , ... , 9 . How- 
ever, they can be resolved. Indeed, if we denote by g a matrix with components 

we can see that due to (1.5) it is Jacoby matrix for Legendre transformation (1.6). 
Matrix (6.8) is non-degenerate since Legendre map A is diffeomorphism due to our 
initial assumptions (see definition 1.1 and ending part of section 1 after it). 
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As we already mentioned above, both (6.6) and (6.7) form a system of homo- 
geneous linear ordinary differential equations with respect to functions r 1 , ... , r" 
and 1 , ... , 9 n . Solutions of these equations form linear space of dimension 2n, 
We denote it by 1. Looking at (6.4), one can see that ip depends linearly on 
r 1 , ... , r n . Hence all time derivatives of ip (i. e. <p, <p, and so forth) depend linearly 
on t 1 , . . . , r™ and on time derivatives of r 1 , ... , t". Due to equations (6.6) and 
(6.7) the latter ones can be expressed linearly through r 1 , ... , r" and 1 , . . . , 8 n . 
This means that for each fixed instance of time t the value of function <p(t) itself 
and values of all time derivatives of this function 

* dt k 

are linear functionals belonging to dual space T*. The dimension of X* is finite, it 
is equal to 2n. Therefore functions <p, <p, <p, . . . , <^( 2n ) are linearly dependent as 
elements of T*. This means that there are some coefficients Co, ... , C^ni which do 
not vanish simultaneously, such that the following equality holds: 

n 

]TC<^=0. (6.9) 

i=0 

For fixed t coefficients Co, ... , Ci n in (6.9) are real numbers. However, if t is 
not fixed, then Co, ... , Ci n depend on t. They also depend on that particular 
trajectory of dynamical system (3.2), for which functions <p^(i) are calculated: 

n 

5^Ci(t) V W=0. (6.10) 

i=0 

Theorem 6.1. For each trajectory p = p(t) of Newtonian dynamical system (3.2) 
all deviation functions ip = ip(t) on this trajectory satisfy the same linear homoge- 
neous ordinary differential equation (6.10) of the order not greater than 2n. 

Saying "all deviation functions", in theorem 6.1 we assume that given trajectory 
p = p(t) can be included into one-parametric family of trajectories by all possible 
ways. Each such inclusion defines some variation vector r = r(t) and corresponding 
deviation function ip = (p(t) on that trajectory. Functions Co, . . . , C2„ in (6.10) 
depend on the trajectory p — p(t), but they do not depend on how this trajectory 
is included into one parametric family of trajectories. 

Theorem 6.1 gives upper estimate for the order of ODE (6.10). For most cases 
this estimate 2n is reached. However, in some special cases real order can be much 
less than 2n. We consider one of such special cases. 

Definition 6.1. We say that Newtonian dynamical system (3.2) satisfies weak 
normality condition if for each its trajectory p = p(t) there is some second order 
homogeneous linear ordinary differential equation 

<p = A(t) <j> + B(t) <p (6.11) 



such that all deviation functions on this trajectory satisfy this differential equation. 
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In other words, weak normality condition is fulfilled when (6.10) reduces to second 
order differential equation (6.11). 

7. Additional normality condition. 

Suppose that weak normality condition is fulfilled somehow. Then all devia- 
tion functions (2.6) arising in shift construction satisfy second order differential 
equation of the form (6.11). Coefficients A and B depend on shift trajectory (see 
definition 6.1), therefore here we should write this equation as 

(p = A(t, y\ . . . , y"- 1 ) if + B(t, y\..., y"" 1 ) <p. (7.1) 

In spite of this minor difference (7.1) is linear homogeneous ODE for ip with respect 
to time variable t. According to definition 2.1, in order to have normal shift we 
should provide vanishing of all deviation functions ifi, . . . , ip n -i in (2.6). Due to 
(7.1) now it is sufficient to provide the following initial data for them: 



<Pi,\ -0, ipA =0. (7.2) 

u=o lt=o 

The only way to provide these initial conditions is to choose proper initial con- 
ditions in (2.1). First part of initial conditions (2.1) says that shift trajectories 
should start the points of initial hypersurface S. We cannot change this condition. 
However, we can specify second part of initial conditions (2.1). Note that at each 
point p € S we have tangent space T p (S) embedded into tangent space T P (M) as 
a subspace of codimension 1. Denote by n = n(p) some nonzero covector that 
vanishes when contracted with all vectors from T p (S): 

n 

(n|r) = ^n s r s = 0, for all reT p (S). (7.3) 

s=l 

It is called normal covector for S. At each point p £ S normal covector n is 
determined uniquely up to a scalar factor. We have no tools for to specify this 
factor canonically, but, nevertheless, we can glue normal covectors from various 
points into a smooth covector-valued function n = n(p) on S (or locally in small 
domains covering all points of S). Now let's compare (7.3) with (2.6) and remember 
that for t = variation vectors T\, ... , t„_i form a base in tangent spaces to S. 
Therefore in order to provide first part of initial conditions (7.1) we should direct 
initial momentum covector p along normal covector of S. This means that we 
should specify initial conditions (2.1) as follows: 



x l \ =x l (p), pA =v(p)-m{p). (7.4) 

lt=0 lt=0 

Here rii(p) are components of normal covector n(p), while v — v{p) is some smooth 
scalar function on S, which is yet undefined. 

When applied to Newtonian dynamical system written in relative form (3.2), 
initial data (7.4) determine initial velocity v implicitly through Legendre trans- 
formation (1.5). However, we can make them explicit initial data if we pass to 
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p-representation in the equations (3.2). Here these equations look like 

x= nWi> Pi = ~n^ +Qi - (7 - 5) 

Equations (7.5) are the equations of Newtonian dynamical system in the form 
relative to modified Hamiltonian dynamical system with Hamilton function H. 
Function H is derived from Lagrange function L as follows: 

" 31 

H = ho\~ 1 , where h = ^v l — - L. (7.6) 

In other words, extended scalar field H is p-representation of extended scalar field 
h, where h is given by formula (7.6). Denominator £1 in (7.5) is p-representation of 
denominator Q in (3.2). It is expressed through Hamilton function H: 

H dpi 

Functions Qi in (7.5) are components of the same extended covector field Q, as 
in (3.2), but transformed to p-representation. Their arguments are x 1 , ... , x n 
and pi, ... , p n . More details concerning formulas (7.6) and (7.7) and Legendre 
transformation in whole can be found in book [2] and in papers [1] , [3] , and [4] . 

Now let's proceed with initial conditions (7.2). First part of these initial con- 
ditions now is fulfilled due to (7.4). We should provide second part of them by 
proper choice of function v = v{p) = v{y 1 , . . . ,y n ~ 1 ) in (7.4). Let's calculate time 
derivative of deviation function ifi using differential equations (7.5): 




s=l s=l 

(7.8 



1 dH 



1 n dx s 

s — 1 s—1 s—1 



In order to calculate time derivatives f? in formula (7.8) we also use differential 
equations (7.5). As a result we get the following formula for f/: 

. s _ d 2 x s d ( 1 8H\ _^dx r d f I dH 
T ' 1 ~ dtdy 1 ~W\tt~dp~s) ~ ^Jhf lte\tt~dp~ s 

+ Qyi [ n dp J - hji Q X T ^ dp J+ Qyi [ft dp J ■ 

Let's substitute this formula into (7.8). Then for time derivative ifi we obtain 

s—1 s—1 s—1 r—1 v / 
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Edp r d ( p s dH\ ^ 1 dH dp s ^ 1 dH s 
W ~dp r \ ^ & dp~s J ~ ^ i ndp~sdy l ~ ^^ds? ' T ' ! ' 

Note that sums in round brackets are identically equal to unity: 

This follows from (7.7). Due to (7.9) two terms in the above formula for ipi do 
vanish. And we get rather simple formula for time derivative tfii : 

El dH dp s 1 dH -J^ 

8=1 fl^^-gfi9? T ' + g ft ^ (7 ' 10) 

If we recall initial conditions (7.2), then from (7.10) we derive 

^1 dH ^1 dH fd Ps \ ^ 



8=1 



4=0 »=1 



Calculating partial derivatives dp s jdy % in the equality (7.11), we should remember 
that p s — v ■ n s . This follows from initial data (7.4). Then 



dps 
dy 1 



dv dn s 1 dv | dn s 

t =o <V ns ^ ^ " cty' Ps v dy*' 



Substituting this expression into (7.11) and using formula (7.7) for 0, we can trans- 
form (7.12) into the partial differential equations for v: 

dv __^v 2 dn s dH " (IdH _ \ s 

dyi- z^n dt d Ps ^ u \ndx° y y Tj ' [< 6) 

Having derived (7.13), we have proved the following lemma. 

Lemma 7.1. Initial conditions (7.2) for deviation functions ipi, ... , ip n -i are 
equivalent to initial data (7.4) for shift trajectories, where scalar function v on 
S satisfies differential equations (7.13). 

In two-dimensional manifold M, when n = 2, hypcrsurfaces are curves. In 
this case we have only one deviation function and only one variable y — y 1 as a 
parameter on S. Therefore (7.13) appears to be ordinary differential equation for 
the function v = u{y). We can set initial value problem 

v{y)\ =*o, (7.14) 

\y=0 

which is always solvable (at least locally) for all vq ^ 0. This means that in two- 
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dimensional case weak normality condition stated in definition 6.1 is sufficient for 
to arrange normal shift of any predefined hypersurface S in M. 

In multidimensional case n ^ 3 situation changes crucially. Here equations 
(7.13) form so called complete system of Pfaff equation for v = v{y 1 1 . . . , y™ -1 ). 
Each separate equation in such system can be treated as ODE. Therefore initial 
condition like (7.14) is the best way for fixing some particular solution: 

i/(po) = vq. (7.15) 

However, initial value problem (7.15) for Pfaff equations (7.13) is not always solv- 
able: some additional conditions should be fulfilled. This is why we had a fork in 
development of theory of metric normal shift in Riemannian geometry (compare 
theses [6] and [7]). This fork is present here as well. 

Definition 7.1. Complete system of Pfaff equations (7.13) is called compatible if 
initial value problem (7.15) for it is locally solvable for all fo ^ 0. 

Let's write Pfaff equations (7.13) formally, denoting by tpi their right hand sides: 

= V i ^y 1 ,...,y"- 1 ) (7.16) 

Due to (7.16) we can calculate mixed partial derivatives of v in two different ways 

Equating (7.17) and (7.18), we get compatibility condition for (7.16): 

0i j {u,y 1 ,...,y n - 1 )=0 ji (v,y 1 ,...,y n - 1 ). (7.19) 

Lemma 7.2. Pfaff equations (7.16) are compatible in the sense of definition 7.1 
if and only if for v ^ left and right hands sides of (7.19) are equal to each other 
identically as functions of n independent variables y , . . . , y n , and v. 

Lemma 7.2 is standard result in the theory of Pfaff equations. Proof of this 
lemma can be found in thesis [6] . 

Definition 7.2. We say that Newtonian dynamical system (3.2) or, equivalcntly, 
dynamical system (7.5) satisfies additional normality condition if Pfaff equations 
(7.13) derived from initial conditions (7.2) are compatible for any hypersurface S 
in M and for any marked point po e S. 

8. Complete and strong normality conditions. 

Both weak and additional normality conditions constitute so called complete 
normality condition in Riemannian geometry (see thesis [6]). We shall keep this 
terminology saying that Newtonian dynamical system satisfies complete normality 
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condition if both conditions from definitions 6.2 and 7.2 are fulfilled 1 . But, apart 
from this complete normality condition, we shall consider strong normality condi- 
tion given by the following definition. 

Definition 8.1. Newtonian dynamical system given by differential equations (3.2) 
or, equivalently, by differential equations (7.5) satisfies strong normality condition 
if for any hypersurface S in M, for any marked point po E S, and for any real 
constant v a ^ there is some smaller part S' of S containing marked point po 
and there is some smooth function v = v{p) in this smaller part S' normalized by 
the condition (7.15) and such that initial data (7.4) with this function v determine 
normal shift of S' in the sense of definition 2.1. 

In simpler words, dynamical systems satisfying strong normality condition arc 
called systems admitting normal shift of hypersurfaces. They are able to move 
normally any predefined hypersurface S in M. 

It is easy to note that complete normality condition is sufficient for Newtonian 
dynamical system to satisfy strong normality condition. Below in section 17 we 
shall prove that it is not only sufficient, but necessary condition as well. 

9. Compatibility equations. 

Now suppose that M is a manifold of dimension n ^ 3. In this case we should 
study compatibility equations (7.19) for Pfaff system (7.13). For this purpose let's 
calculate partial derivatives (7.17) and (7.18) explicitly: 

d 2 v v 2 dH d 2 n s (LUK-O^ d2x ' S 



dy 1 dy 3 j-— J Q dp s dy l dy 3 J— J \ VL dx s J dy 1 dy 3 

2^2^ 02 f)n. f)n„ fh.i f)„3 +l/ [ci < fir* Qs O fl~r Qr T * T 3 4 



—J ^-J \ £1 2 dp s dp r dy 1 dy 3 \ ft dx s 8 J \ ft dx r 

2 v 2 ( 1 dH n \ dH dn r „ , v 2 dH ( 1 dH „\ dn s r 



+ n \n dx s Q j dp r dy 3 Tl + n d Ps \n dx r Qr ) d y i Tj 
hihC d P* \ n d vr) dy 1 9y 3 ~ hihi" dps \ n dxr / dyt T] 



2^/1 ®H \ ® n r s 9/1 dH 

s=l r =l ^ fi/y s=1 r=1 \ / 

For Pfaff equations (7.13) to be compatible, right hand side of the above equality 
should be symmetric in indices i and j. Some terms there are obviously symmetric. 
Below we shall not write such terms explicitly denoting them by dots: 



d 2 v 



vvV— (-^L-o\ — - 2 9 ( l dH X\ drir s 

~' + i Qs ) dp~r V ~dx* \U ~dp~ r J J 'dy 3 T ' 1 + 



dyi dy 3 , , 



1 In two-dimensional case n = 2 complete normality condition reduces to weak normality con- 
dition since additional normality condition in this case is always fulfilled. 
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" " f^dHfl dH_ _ \_ 2_d_f}_dl£_ n \\dn 1 r 

2^2^ { Q dps [q dx r V ^ VrJ J Qyl Tj 

8 fl dH\ dv dn r d f 1 dH \ dv 



s—1 r—1 



§ 5 vPs dp s V n dPr ) dy l dyi ^^ Ps d Ps \n dx r Qr ) d v i T] 

S — X 7" 1 — X S — X IT — X 

v^v^ 3 <9 / 1 dH \ dn s dn r d fl dH \ s r 

S — X 7* — X S — X IT — X 

In the above calculations we used formula (7.12) in order to express partial deriva- 
tives dps/dy 1 through dn s /dy l . As a result we have got partial derivatives dv/dy l 
in the above expression. In order to eliminate them now we shall use (7.13): 

d 2 v_ fvf-(i--ol-- 2 9 (- dH \ 



dy % dyi 7^~i7~~i \ ^ ® xS ) ®Pr dx s \ Q dp, 

^ 2 8(1 8H\/1 8H \\ dn r s / v 2 dH 



m \ 2^_(}_dH__ n \ v 2 p q d fl dH ] 

ndx- Qr ) v dpX^idx- Qr ) + 2^ n d Pq \ndxr Qr l x 



aff\an, r " " / " v A Vq d f 1 dH\ dH v A dQ dH 

x dp~ s ) W T3 \ ^i~^~dp^\Jidp^) Ws + wdp~ s dp~r / 



X 



d^dnr yy/A d fl dH \fl8H 

WW + ^ 1 ^^ Pq WAn^~ Qr )\ndx^~ 



'» 3 



Term with product of partial derivatives dn s /dy l and dn r jdy^ in the above expres- 
sion can be replaced by dots. Indeed, one can easily check up that coefficient of 
such product of derivatives is symmetric in indices s and r: 

A v 3 Pq d f 1 dH\ dH v 3 dn dH <A v 3 p q dn dH dH 
£1 dp q \Q dp r J dp s Q 2 dp s dp r ^--J Q 3 dp q dp r dp s 



v'p q d 2 H dH V "Pq 
fl 2 dp q dp r dp s ^ fl 2 dp q dp s dp r 1 fl 2 dp s dp r ' 



Ev p q urn un x ^ 
(( =1 ~W dp q dp r d P ~s + q ^i 



Now let's study the term with t? tJ. For the coefficient in this term we derive 



^ a fi dH \f i dH \ d fi dH 



3=1 
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~ hi ^ ^ W q dx* : + n d Pq dx r ~l^)\nd^~ Q ' 



1 dn dH t 1 d 2 H dQ r \ _ , j^vpq d 2 H dH 



CI 2 dx s dx r Q dx s dx r dx s I ' ^-f fl 2 dp a dx r dx s 

/ q=l 111 

E vpq dQ r dH y> vp q dH dH Y^y^ ^PgPfc d 2 H dH 
~fT~dpZdx~ s ^lr r & P ~lhf Qs : + 2^2^ ^2 d d dx^ Qs ^ 

q=l H q=l H q=l k=l H 

Evp q d 2 H ^ dQ r ^ vp q d 2 H dH v d 2 H 

IT d Pq dx r Qs + 2-; ypq ~dpl Qs + 2^ IF d Pq dx s d&~Q. dx r dx s + 

q=l H 5=1 H 5=1 4 

dQr _ ^ vpq dQr dH ^ vp q dH dH v VgVk 

dx s 2-<n d Pq dx s 2^ n 2 d Pq dx r Qs + 2^2^ fi2 x 

5=1 ^ 9 5=1 ^ 9 5=lfe=l 

d 2 H dH "vp q d 2 H + y vp ^Rlq +v ^9l 

dp q dpk dx r s 2-^ ft dp q dx r s ~^^ q dp q s V dx s ' 

In the above calculations we used formula (7.7) for Q,. As a result for partial 
derivative (7.17) we have derived the following expression: 



Z^Z-^ l O \ O ft^s ^ s f)r>„ \ O f)r>„ I Z-^ 



Qyi dyj Z-j Z_j y fi Q x s J dp r dx s \ Q, dp., , 



X 



d (I dH\(l dH ^\\ dn r s s^s^(v 2 dH ( 1 dH ^\ 

x Pq d Pq {ndpr)\n dx* Qs ) ) d y > T% + ^ ^ \n d Ps [n dx- Qr ) 

a d (\ dH \ ^i^pq _d_f }_dH_ _ n \ dH_\ dr^ r 

v d Ps \n dx r Qr ) + 2^ n d Pq \n dx r Qr ) d Ps ) dt Tj + 

Y" V | y v P<i dC * r dH \ " ^ ^ ^ g n \ \^ 

+ 2^ 2^ I 2^ n a?3 + Z-j ^2 9 dx r ^' + 2^ 2^ 

s=l r=l \ 5=1 ^ q 5=1 ^ g 5=1 fe=l 



fi 2 

5=1 * 1 5=1 * * 5=1 fc=l 

J^HdH ^ V p q d 2 H ™ 9Q 9Q r \ s r 

Next step in transforming the above expression is based on the following equality: 

s=l r =l y s=l r=l y s=l r=l y 

Here, as we already used above, we denoted by dots terms symmetric in indices i 
and j. Further for our particular ct rs and Pr S we derive 

r Rr v 2 (\dH \dH 2 d/ldH\" 2 

a °-^ = n\sid^- Qs )Wr- v dx-\ndp-) + ^ PqX 
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9p 9 9p r y ^ft dx s Qs ) d Pr \n dx s Qs ) 4^ 



d fl dH Q\dH_y^y^v 2 pp d 2 H 8H 
dp q \ £1 dx s 7 9p r 4rf r~f ^ 2 P<? Pfc ^P<? ^Pfc dPr 



4^ ^ ^ 9 <9p<j dp r s V dp r ^ ^ 9 ^Pg ' 



Summarizing all above calculations, for partial derivatives (7.17) we obtain 



d 2 v 

Qyl Qyj 

n 9 
9=1 



n n o 



d 2 g n 2 9Q S 
a — Qs-v 77— 

Op q Op r Op r 



Ev 2 dH dQ s \ dn r . 
77 — — I — t/ + 



9=1 



f2 <9p r 9p 9 / dy-? 



+ z^/M A, n &e* ^ ft 2 d Pq dx r ^ s ^ L> ' 



9=1 



d 2 H dH 



dp q dpk dx 



un_ „ _ ^ VPq d 2 H 



^P? Qs+ V 



9=1 fe=l 



f! 2 



9=1 



dPq 



dx s 



T- T ■ . 

* J 



(9.1) 



In a similar way for partial derivative (7.18) one can get analogous expression: 



n n / n n 9 



dyi dy l 



+ EE EE 

s=l r=l \q=l k=l 



d 2 H dH 



Qr 



Ev 2 d 2 H 2 dQ r ^ v 2 dH dQ r \ dn s r 

__ i n Pq dp~ q ~dp~ s Qr ~ v ^ + ^ l n Pq dVs~d P ~)W Ti + 

y^\p / \ - vp q dQs oH \ ^ dH ^PgPfc 

, 1 I 1 " 9 P9 ^ , fi2 9 ?9 feS 1 , & 

s=l r=l \ q=l rH q=l rH q=\ k=\ 

x J^_dH_ Qr _y^P 1 ^H_ Qr + y iyp dQ s | ^ ggA _ s 



(9.2) 



Now we are able to write compatibility equations (7.19) in explicit form. It is 
sufficient to equate partial derivatives (9.1) and (9.2) to each other. However, 
we shall not do this right now because we need some additional information in 
order to treat arising compatibility equation properly We should replace partial 
derivatives in (9.1) and (9.2) by covariant derivatives, and we should understand 
some facts from theory of hypersurfaces in a manifold equipped with the only 
geometric structure given by Lagrange function L. 
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10. Differentiation of extended tensor fields. 

Let's consider smooth extended tensor fields in their p-rcprcscntation as given 
by definition 5.1. They form a graded ring with respect to standard operations of 
addition and tensor product. We denote it as follows: 

oo oo 

T(M) = 007J(M). (10.1) 

Graded ring (10.1) is equipped with operation of contraction, which is also standard. 
Moreover, (10.1) possesses the structure of algebra over the ring of smooth functions 
in T*M. For this reason it is called extended algebra of tensor fields. 

In extended algebra of tensor fields (10.1) one can define canonical covariant 
differentiation, which is called vertical gradient or momentum gradient: 

V: TJ(M) — ► TJ +1 (M). (10.2) 

In local chart momentum gradient (10.2) is determined by the following formula: 

QX^ " ir 

V'.X]' [ — (10.3) 

^■■■Jo dp q 

To define another covariant differentiation in T(M) one need some additional geo- 
metric structure in M. This is so called extended affine connection. We shall not 
discuss the nature of this structure (see thesis [6] ) . Note only that in local chart it 
is given by its components 1*-, which obey standard transformation rule: 

n n n n f) r F' m 

r "j = E E E sk m T * T I f "c + E^^r- ( 10 - 4 ) 

m—la—lc—l m— 1 

Here 5* and Tj are components of transition matrices S and T binding coordinates 
x 1 , . . . , x n and x 1 , ... , x n in two overlapping local charts of M: 

S) = T\ = (10.5) 

3 d& 1 dxi v ' 

Unlike components of standard affine connection, components of extended affine 
connection V in p-representation depend not only on x 1 , ... , x n , but also on com- 
ponents pi, . . . , p n of momentum covector p: 

It !,/•<•' '•"•/'! Pn)- (10.6) 

Thus, if M possesses extended affine connection T, one can define horizontal gra- 
dient or spatial gradient V in T(M). In local chart it is expressed by formula 



dx n-ir JLJL dX\ 



yj y»i-V _ 31- 3s \ S^S^„ I '•• ""31-3 
V\, -3 ~ + 2s l^ PaL lb' 



Q - lfc=1 (10.7) 

r n s n V ' 

E\ ^ pifc y-ii... a k ... i r _\^\ 1 pfcfc vil ir 
9 a fc ji 3s <7Jfc ji- 6*- J« " 

fc=lo fc = l fc=16 fc = l 
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Horizontal gradient V increments by 1 the number of lower indices of tensor field: 

V: T:(M)^T: +1 (M). (10.8) 

Using Legendre map, one can transform extended tensor fields to v-representa- 
tion. In v-representation vertical gradient V is called velocity gradient. It is a map 

V: t:(M)^T: +1 (M) (10.9) 

similar to (10.8). In local chart this map (10.9) is given by formula 

ox i.-i f 

Two versions of vertical gradient (10.3) and (10.10) can be related to each other by 
means of formula using components of matrix (6.8): 

n 

% ( x l:X A > X g gk v*xj; ;;; ; : o a. (io.il) 
fc=i 

Remember that inverse Legendre map in local chart is given by explicit formula 

v* = ^ (10.12) 

Therefore we can introduce matrix g _1 with the following components: 

8 2 H 

o" = a^i (iai3) 

Matrices (6.8) and (10.13) are inverse to each other when taken in the same rep- 
resentation, i. e. when both brought to p-representation or when both brought to 
v-representation. Matrix (10.13) is used in formula, which is similar to (10.11): 

n 

V"(*ri o A- 1 ) = ]T V k Xt:}: . A" 1 . (10.14) 

fe=i 

Looking at formulas (10.11) and (10.14), we see that matrices (6.8) and (10.13) 
here do part of work that metric tensor does in Ricmannian geometry 

In order to define horizontal gradient for extended tensor fields in v-representa- 
tion we should transform connection components (10.6) by means of Legendre map: 

r k j ^T k j °\ = T k j (x 1 ,...,x n ,v 1 ,...,v n ). (10.15) 

Then, using (10.15), we can define horizontal gradient by a formula in local chart: 

r\yii...i r n n ayii-..i r 
^VSl-Js Q x q Z^2^ V l W dv b 

a=lb=1 (10.16) 

r n s n V ' 

I \ " \ 1 pifc v-i\...a k ... i r _ \ y \ y -p6 fc yil ir 

' 1 '\/ is 2^ Qikii-b k ...i 3 - 

k=l a k = l fe=16 fc =l 
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In spite of transformation rule (10.15) for connection components, in general, (10.7) 
and (10.16) are not different representations of the same tensor field. We have the 
following equality binding these two horizontal gradients: 

n 

\U.y:.; ;; ; A} ^.v: :;; ; a ^y; ;:: ; a. (10.17) 

However, there is special case, when (10.7) and (10.16) do coincide. This is when 

%V S L = 0. (10.18) 



In this case (10.17) turns to equality, which is similar to the equality (10.11): 

%(Xj-;j : oA) =V q X£;;X°\. (10.19) 

Definition 10.1. Extended connection T is called concordant with Lagrange func- 
tion L if the equality (10.18) holds. 

For concordant connections due to (10.19) horizontal gradient can be calculated 
either in p or in v-representation yielding the same result, but in different variables. 
Note that the equality (10.18) can be replaced by equivalent equality for Hamilton 
function H in p-representation. It looks like 

W q V s H = 0. (10.20) 

This equality (10.20) can be derived from (10.18) by direct calculations. 



11. Differentiation along lines and hypersurfaces. 



Let p — p(t) be a smooth parametric curve in our manifold M. In local chart 
with coordinates x 1 , ... , x n it is represented by functions 



' x 1 ^x^t), 
c n = x n (t). 



(11.1) 



Suppose that at each point p of this curve (6.1) some tensor of the type (r, s) is given, 
i.e. we have tensor- valued function X = X(i). In local chart this tensor function 
is expressed by its components X 1 ^"^ = A^'"^(i). If we had standard affine 
connection in M, we could define covariant derivative of tensor function X(t) with 
respect to parameter t along the curve. This is another tensor function Y = VfX 



with components Y- 



1 ji- 



given by the following formula: 



dx l:Z 

dt 



r n n 



k ... i r 

is 



k=l m=l a k = l 
s n n 



(11.2) 



EV \ ^ ™m -p bk x il ir 
mj k yl j 1 ...b k ...j s - 

k=l m=l b k = l 
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However, with extended connection T we need come additional data in order to 
apply formula (11.2) to tensor function X = X(t). Indeed, looking at (10.6) and 
(10.15), we see that in p-representation we need functions 

pi(t), ... , Pn (t), (11.3) 

while in v-representation we need another set of functions: 

v\t), ... , v n {t). (11.4) 

Taken by themselves, functions (11.3) and (11.4) arc components of covector- 
function p = p(i) and vector-function v = v(t). But combined with functions 

(11.1) , they define lift of initial curve p = p(t) from M to cotangent bundle T*M 
and to tangent bundle TM respectively. Once such lift is given, we can use formula 

(11.2) . Covariant derivative given by this formula then is called covariant derivative 
with respect to parameter t along the curve p = p(t) due to its lift q = q(t). 

Now let's discuss the question of how tensor function X = X(t) on curve could 
be defined. Surely, it can be given immediately as it is. However, often tensor 
function X = X(f) comes to curve from outer space, i. e. from manifold M. For 
example, if in M or at least in some neighborhood of our curve some standard 
(not extended) tensor field X is given, we can restrict it to the curve p = p(t) thus 
obtaining tensor function X = X(i). In the case of extended tensor field X we 
cannot restrict it to the curve p = p(t) immediately. We should first choose a lift of 
this curve q = q(t) in TM or in T*M , then we could restrict X to the lifted curve. 
Now suppose that tensor function X = X(£) on the curve p = p(t) is obtained from 
extended tensor field X in this way. Then covariant derivative (11.2) for X(t) is 
determined. By direct calculations we can prove that this covariant derivative can 
be expressed through covariant derivatives (10.7) and (10.3) in p-representation: 

n n 

v^:::}: = £ v k x>-± ■ ± k + £ v**£::£ ■ v tPk . (11.5) 

fe=l k=l 

In v-representation we have similar expression 

n n 

^Xf.l = £ V^;;;} : • i k + £ V^;; j; • V,^ (11.6) 

k=l k=l 

that uses covariant derivatives (10.16) and (10.10). Here V t v k and \7 t pk are com- 
ponents of covariant derivatives for vectorial and covectorial functions with com- 
ponents (11.4) and (11.3) that determine lift of curve p = p(t) to TM and T*M 
respectively. From general point of view, formulas (11.5) and (11.6) are nothing, 
but the rule of differentiating composite functions written in terms of covariant 
derivatives. They are obtained by direct calculations. 

Now let S be a smooth hypersurface in M. Denote by y 1 , . . . , y 11 ^ 1 coordinates 
of point in some local chart of S. Then the following smooth functions 

' x 1 = x 1 (y 1 , . . . ,y n ^ 1 ), 

< (11-7) 

k x n = x n (y 1 ,...,y n - 1 ) 
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represent hypersurface S in local chart of M. Partial derivatives of functions (11.7) 
determine tangent vectors n, ... , r n _i to S as it is done in (2.5). Taking one vari- 
able y % in (11.7) and fixing all others, we can consider (11.7) as n— 1 families of func- 
tions of one variable, each corresponding to some 
family of curves on S. These are coordinate curves 
forming coordinate network on S (see Fig. 11.1). 
One can lift coordinate curves to cotangent bundle 
T*M by means of normal covector n = n(p). How- 
ever, we shall use another lift given by momentum 
covector (see initial data (7.4)): 

p = u(p) ■ n(p). (11.8) 

Then applying Legendre transformation (1.5) to 
covector (11.8), we obtain vector 

v = v(p)=VI (11.9) 

at each point p G S. Vector v(p) is transversal to S 
as shown on Fig. 11.1. This follows from regularity 
of Lagrange function (see definition 1.1). Now, using lifts defined by (11.8) and 
(11.9), we can apply formula (11.2) to coordinate lines of hypersurface S. 

Remember that t = y l is a parameter of i-th coordinate line, while t% is tangent 
vector to this line corresponding to parameter y l . Covariant derivative V t with 
respect to parameter t = y % along i-th coordinate line by tradition is denoted by 
V Ti . This derivative can be applied to any smooth tensorial function defined at 
the points of hypersurface S. Let's apply it to vector- function Tj. As a result ge 
get another vector-function V Xi Tj on S. Similarly one can consider vector- function 
V Xi v and covector-function V Ti p on S. The latter one appears to be most important 
for us. It will be studied in section 13 below. 

12. Covariant form of compatibility equations. 

In section 7 we have derived Pfaff equations (7.13) for scalar function v on 
hypersurface S. Then in section 9 we studied compatibility condition for these 
Pfaff equations. Explicit form of compatibility condition could be obtained by 
equating partial derivatives (9.1) and (9.2). But as a result we would obtain huge 
equality difficult to observe. This means that formulas (9.1) and (9.2) require some 
preliminary transformations in order to simplify further analysis of compatibility 
condition they lead to. 

First of all we are going to replace partial derivatives dn T jdy^ and dn s /dy l 
by covariant derivatives V Tj p r and V Ti p s , assuming that we have some symmetric 
extended affine connection in M. Let's turn back to the equalities (7.12) and (7.13), 
let's substitute dvjdy % from (7.13) into (7.12) and remember formula (10.12): 



d P g V" fxr v r p s \ dn r ^ / 1 dH \ r 



(12.1) 
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Denote by PJ coefficient in front of partial derivative dn r /dy % in formula (12.1): 

pr = 5l-^. (12.2) 
Formula (12.2) for PJ can be rewritten in two equivalent forms: 

PI = s:- p -^2L = s:-q^ (i2.3) 

(compare with formulas (10.12) and (1.5) above). Now we see that (12.3) are compo- 
nents of operator- valued extended tensor field P either in p and in v- representations 
respectively. It is easy to check up the following identity: 

P 2 = PoP = P. (12.4) 

Formula (12.4) means that P is a projector- valued extended operator field in M. It 
projects along velocity vector v onto a hyperplane defined by momentum covector 
p. When restricted to hypersurface S due to lift (11.8) in p-representation or due to 
lift (11.9) in v-representation, it becomes projector field that projects onto tangent 
hyperplane to S along velocity vector v. In this restricted form components of P 
appear in formula (12.1). This formula now can be written as 

Let's use the following identity that can be checked up immediately: 

E^. = E(*S-^)p. = 0. (12.6) 
Applying (12.4) and (12.6) to the equality (12.5), we find that 

Vp : ^ = 1Vp : |. (12 . 7) 

^ s dy l v ^ s dy l y J 

r=l y r=l y 

Looking at (12.5), we see that we cannot express partial derivative dn s /dy l in 
a pure form. However, as appears, formula (12.7) is sufficient for our purposes. 
Indeed, the term containing partial derivative drir/dyi in (9.1) can be written as: 



n n n n 



Therefore we can use the above identity (12.7) in order to transform this term: 
T^i 7~i £i f~i V n d Pi d Pk d Pq J dyi 
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As a result of this transformation for the whole expression (9.1) we obtain 

d 2 v A\A\ /AA v_ d 2 H dH 

dyi dyj - ■ ■ ■ + L, I, [L, L, 02 ft m ^ apr «• 

Ev_ d 2 H _ dQ s t A\ v dH dQ s \ dp r 
n P9 a P9 a Pr ^ s 17 s Pr + ^ si Pq d Pr d Pq j d y > Tl + 

(12.8) 



n n 



s =l r =l \ q=l '1 q=l 1 1 q=lk=l 

d 2 H dH 



n _ \A vj>q_ d 2 H A\ ag^ agA . , 
9 Pg a Pfe ar- Vs ^ n a Pg dor y * + Z^ft dPq ^ s + v dx° ) T * T > 



Now we can replace partial derivative dp r /dy j in (12.8) by covariant derivative 
V Tj p r - For this purpose we shall use formula (11.2) with t = yK This yields 

rj n n 

y a =lo-=l 

Let's substitute (12.9) into (12.8) and let's use formulas (10.3) and (10.7) in order 
to express partial derivatives of H and Q s through covariant derivatives: 

~9 n n n 

a v 



s=l r =l o=l ^ ' 



y a s=l r=l q- 

+ EE n2 ^~Eft^ Qs-^ r Q s + (12.10) 

s— 1 r— 1 \ q—1 

+ Eft ^ + v -^f- Qs + E ^ft ^ ^) < r l- 

9=1 "" 9=1 / 

By transforming (9.2) we can obtain another expression for the same derivative: 

-9 n n n 

d v - — - — — 



dyi dy % 



■■•-EEE(p'| + ^q) v PI V TiPs rj + 

g=l r=l q=l ^ ' 



+ E E i 1 ^^- Qr^q Qr - V sQr + (12.11) 

s— 1 r— 1 \ " q—1 

A 1 ^ v,i? - v v,h A\ - \ 

+ Eft^^ V "^ + — jf-Qr + 5>P g QrV*Q a r/rj. 

9=1 9=1 / 

Here, as in section 9, we denote by dots terms symmetric in indices r and s. They 
do not affect ultimate compatibility equation, which will be derived from (12.10) 



2d 



RUSLAN A. SHARIPOV 



and (12.11). In the equalities (12.10) and (12.11) we used matrix = H 
from (10.13) as metric tensor and introduced the following notations: 

n n n 

p q = Y.9 qk p<n \p\ 2 = J2Y,9 qk p«p»- ( 12 - 12 ) 

k=l q=l k=l 

In general, matrix (10.13) is not positive, therefore |p| 2 in (12.12) is not necessarily 
positive number for p ^ 0. 

Let's equate right hand sides of (12.10) and (12.11). As a result we obtain 
compatibility equation (7.19) written in terms of covariant derivatives: 

n n n 



n n n / ^ \ 

EEE + ^Q s ) P r q S7 TjPr rt 

s = l r=l q=l ^ 



n n n 



s— 1 r— 1 q— 1 



1^ 1^ — 02 Qs ~ l^P* o Qs ~ VrQs + 



fi 2 ^° n 

s— 1 r—1 \ q—1 



(12.13) 



+ X> ^ 0s + T Qa + ^ ^ s ^ j r/ rj- 

9=1 "' 9=1 / 

[ f]2 Qr-2_^P<1 Q Qr-y a Qr + 

s—1 r—1 \ q—1 

+ ltpi^f + ^Qr + J2Pl ®r V 9 Q a ) 7? r] . 

9=1 9=1 / 

The most important point here is that the form of equation (12.13) does not depend 
on which particular connection is used in covariant derivatives. One can choose any 
symmetric extended connection r in M. 

13. Geometry of hypersurfaces. 

Though most terms in compatibility equation (12.13) can be treated as compo- 
nents of extended tensor fields in M, the equation in whole is related to hypersur- 
face S. This reveals if we look at partial derivatives V Tj p r and ^ Ti Ps- For further 
treatment of (12.13) we should learn how to treat these partial derivatives. 

Let r be some arbitrary tangent vector to 5. It represented as linear combination 
of basic tangent vectors: r = a 1 • n + . . . + a™ -1 • r„_i. Let 

n n— 1 

f(r) = V T p = E E ( aJ V ^Pr) ■ dxT - I 13 - 1 ) 

r=l j=l 

Then (13.1) defines linear map /: T P (S) — > T*(M). Let's consider composite map 

b = -P*o/oP. (13.2) 
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Projection operator P* in (13.2) is a conjugate operator for projector P with com- 
ponents (12.2). Remember that P projects onto the subspace T p (S) in T p (M) for 
p G S. Therefore linear map b : T p (M) — > T*(M) is correctly defined by formula 
(13.2). Each map from T p (M) to conjugate space T* (S) is given by some bilinear 
form. In our case this bilinear form is defined by formula 

6(X,Y) = (b(Y)|X>. (13.3) 

Due to the presence of projection operators P and P* in (13.2) we have 

6(X,Y) = &(P(X),Y) = &(X,P(Y)). (13.4) 

Theorem 13.1. Bilinear form (13.3) defined by linear map (13.2) is symmetric. 

Bilinear form (13.3) restricted to tangent space T p (S) of hypersurface S is known 
as second fundamental form of hypersurface S. Its components 

/3 ij =b(r i ,T j ) (13.5) 

define tensor field in inner geometry of hypersurface 5*. But for our purposes coor- 
dinate representation of bilinear form (13.3) in outer geometry is more preferable: 

n n 

b = V" ^2 hj dx % dx j . 

i=l o = l 

In order to prove theorem 13.1 we need some preliminary results. First is given by 
the following lemma for extended connection components. 

Lemma 13.1. For any symmetric extended connection T in M and for any fixed 
point q = (pcbP) of cotangent bundle T*M there is local chart of M such that 
all connection components T^(q) = T^x 1 , . . . , x n ,p\, . . . ,p n ) in this chart become 
zero at the point q = q . 

Lemma 13.1 is formulated for extended connection in p- representation and for 
its components (10.6). Similar lemma can be stated for extended connection in 
v-representation and for its components (10.15). 

Lemma 13.2. For any symmetric extended connection T in M and for any fixed 
point q = (po, v ) of tangent bundle TM there is local chart of M such that all 
connection components T^-{q) = T^x 1 ,... , a;™,?; 1 ,... , v n ) in this chart become 
zero at the point q = q . 

Here one should emphasize that lemmas 13.1 and 13.2 assert vanishing of T^- 
only at single point q = q . Sometimes they can vanish in whole neighborhood 
of this point, but this is not the case of general position. Proof of lemmas 13.1 
and 13.2 is rather standard. It is based upon formula (10.4). One can find proof 
of analogous lemma in Chapter V of thesis [6]. With minor changes this proof is 
applicable to lemmas 13.1 and 13.2. 

Now suppose that x , ... , x n are coordinates in a chart, existence of which is 
asserted by lemma 13.1. Then, using transformation formula (10.4), one can check 
up that linear change variables (i.e. that, for which transition matrices S and T 
in (10.5) are constant) preserves the property of vanishing of r* 'Aq) at the point 
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Qo = (pojP)- If we apply this fact to a fixed point po £ S, we can choose local 
coordinates x 1 , ... , x n in some neighborhood of po such that 

1) fixed point po is the origin in coordinates x , . . . , x n , i.e. x l (po) = for the 
whole range of index i — 1, ... , n; 

2) hypersurface 5 is given parametrically by the equations 

x 1 = y n , 



x n-l = y n-l } (13.6) 

ar" = 0(y 1 ,... ,y"- 1 ), 

where ^(y 1 ,... , J/™" 1 ) is some smooth function of parameters y , ... , y™ _1 
vanishing at the origin, i. e. z(0, ... , 0) = 0, and having extremum there; 

3) normal covector n = n(y 1 , . . . , y" -1 ) of hypersurface S related with momentum 
covector p on S by the equality (11.8) is given by its components: 

dz 

n=(-z[, ...,-4,1), where z[ = — ; (13.7) 

4) connection components = r^(p,p(p)), where p = p(p) is determined by 
the equality (11.8), do vanish at the origin p = pq. 

Formula (13.7) is derived directly from (13.6). Indeed, if we calculate tangent 
vectors n, ... , r„_i by using formula (2.5), for t£ we obtain 

for s^i, n, 

1 for s = i, (13.8) 
z- for i — n. 

Now it's easy to see that covector (13.7) is orthogonal to vectors n, ... , r„_i in 
the sense of the equality (7.3). As we know, normal covector of hypersurface is 
determined up to a scalar factor. This uncertainty in (13.6) is eliminated by the 
condition that last component of n is equal to unity. 

Let's consider components of covector V x p in (13.1) for those special coordinates 
x 1 , . . . , x n we have chosen above. Applying (11.2) and (12.7), we derive 

n n 



Applying formula (11.2) again, for covariant derivative Vin r we obtain 

dn 
dy 



V r = ^-EECn a r/. (13.10) 

a— 1 a— 1 



Remember that in the above formula T^ r = for p = p . Therefore, taking into 
account formula (13.7), from (13.10) we derive the following equality: 

f z'f for r < n, 
V Ti n r \ = r (13.11) 
\p= Po (0 for r = n. 
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Projector P projects onto the hyperplane T P (S) and Tj € T p (S), hence P(t^) = Tj. 
Then for b(rj), applying (13.1), (13.2), and (13.9), we derive 

n n 

bfa)| = -£E>^ V Ts n r )-dx s . (13.12) 
Further, using formulas (13.3), (13.5), (13.11), and (13.12), wc obtain 

n n n n 

&A_ =-EE< v ^ T i=-EE" p » z "r T i' ( 13 - 13 ) 

Now remember formula (12.2) for PJ and formula (13.7) for n. Moreover, let's 
remember that function z(y 1 , . . . ,y n ^ 1 ) has extremum at the origin. This means 
that z[ = for p = p . Then from (12.2), (13.7), and (13.8) we get 

v r 5 n 

p s\ = 6 s-^r> T °\ = d °- ( 13 - 14 ) 

lp=Po u lp=Po 

Substituting (13.14) into (13.13) we now obtain the following equality: 

H=„ = = w <13 - 15) 

Looking at (13.15), it's easy to see that second fundamental form of hypersurfacc 
S is symmetric: fyj — [3ij. This is in concordance with classical results for hyper- 
surfaces in Riemannian manifolds. But now we are in quite different geometry. 

Second fundamental form (3 with components given by formula (13.5) is a tensor 
field in S. It is obtained by restricting quadratic form (13.3) to inner geometry of 
S. Let's calculate components of this form in outer geometry, i. e. in local chart 
of manifold M at the point p — po- Due to (13.14) tangent vectors n, . . . , t„_i 
coincide with n — 1 coordinate tangent vectors at the point p : 

d 

tA =E i = — for i = l, ... ,n-l. (13.16) 

\p=pa UJj 

As a consequence of (13.16), (13.3), and (13.5) we obtain 

bn = b(E i ,E j ) = b(n,Tj) = fc 3 for 1 sC i, j < n. (13.17) 
Let's take n-th coordinate vector E„ = d/dx n . From (13.14) we derive 

n—li n—1 



p(E„)i =E^- E ^ = -Ef - E - = -Ef -Tr- 

lp=P0 i=1 r=1 V n i=1 V„ 

Using the above equality and formulas (13.4) and (13.5), for i < n we find 

b in = b(E u E n ) = b(Ei, P(E„)) = - E (13.18) 



r=l 
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In a similar way for component b ni of quadratic form (13.3) we obtain 

b m = 6(E n , Ei) = 6(P(E n ), EO = - 2 ( 13 - 19 ) 

r— 1 

Comparing (13.18) with (13.19) and taking into account symmetry of (see 
(13.15) above), we find that bi n — b ni for i < n. Combining this equality with 
(13.17) and taking into account symmetry of again, we get 

bij = bji (13.20) 

for all i and j. We have proved the equality (13.20) at one point p = po by choosing 
special local chart in M. But we can choose arbitrary point of M for po- Besides, 
bij are components of tensor. Therefore, being symmetric in one local chart, they 
keep symmetry when transformed to another chart. This proves theorem 13.1. 

Theorem 13.2. Let q = (po,p) be some fixed point of cotangent bundle T*M 
with p ^ and let projector P be the value of projector-valued extended tensor field 
(12.4) at this point. Then any symmetric quadratic form b in T po (M) satisfying the 
equality (13.4) can be determined by some hypersurface S passing through the point 
Po and tangent to null-space of covector p at this point. 

Proof of theorem 13.2 now is very simple. Indeed, quadratic form b satisfying the 
equality (13.4) is completely determined by its restriction to null-space of covector 
p. Then we can choose local chart in M with r^(p ,p) = and with first n—1 
coordinate vectors Ei, ... , E„_i all being in null-space of covector p at the point 
po. Let's define matrix (3 by means of components of quadratic form b: 

[3 l3 = b l3 = b(Ei,-Ej)\ for 1 < i,j < n. (13.21) 

lp=Po 

Now it's sufficient to define hypersurface S by parametric equations (13.6) and take 
the following function z — z(y 1 , . . . , y n ^ 1 ) in them: 

n n 

z = -j- EE&^V- (13.22) 

»=1 3 = 1 

Here u = v{po) is a constant from (7.15). Due to equalities (13.15), (13.17), and 
(13.21) it is clear that such hypersurface S reproduces quadratic form b used to 
define it through (13.21) and (13.22). Thus, theorem 13.2 is proved. 

14. Additional normality equations. 

Now let's apply theorem 13.2 to the study of the equation (12.13). Let's fix some 
point po in M and some covector p at this point. Then, relying upon theorem 13.2, 
let's take hypersurface S passing through this point tangent to null-space of covector 
p and such that its second fundamental form (3 is zero at the point p = p . This 
means that quadratic form b in (13.3) is also zero for p = p . Then due to (13.1) 
and (13.2) we get V Tj p r = 0. Equivalently, V Ti p s = 0. Substituting these two 
equalities into (12.13), we find that first two sums in (12.13) do vanish. Looking at 
other terms, we see that they are components of an extended tensor field of type 
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(0, 2) contracted with two vectors T\ and Tj tangent to S. These two vectors depend 
on the choice of parameters y 1 , . . . , y n , i. e. on the choice of local chart on S. By 
choosing this local chart properly we can associate Tj and Tj with two arbitrary 
vectors in T Po (S). This means that we can write 



Ti=P(X), tj =P(Y), (14.1) 

where X and Y are two arbitrary vectors in T po (M) . Due to arbitrariness of vectors 
X and Y in (14.1) from rest part of (12.13) we derive 

EE fl2 Q'-Y.P*— — Qs-VrQs + 

8=1 r=l \ 3=1 



+ I> ^rr + 7T^ + Ep^ s V*Q r ] p l = 



3=1 3=1 



(14.2) 



^^f |p| 2 v s h ^ " \7 s vm 

= E E I Q2 Qr - 2^Pq n Qr - VsQr + 

s—1 r—1 \ q—1 

+ X> ^ qQr + "7T 0r + E>* Qr ^ qQs ) p/ p ^ 

Note that the equations (14.2) are partial differential equations for components of 
extended covector field Q, which is used in the equations of Newtonian dynamics, 
when they are written in the form (7.5) relative to modified Hamiltonian dynamical. 
They are written in terms of covariant derivatives (10.3) and (10.7). Though we 
used some special hypersurface S in order to derive them, in their ultimate form 
they do not depend on any particular choice of S and, moreover, they can be written 
in the absence of S at all. 

Now remember that P(rj) = t% and P(tj) = Tj. Therefore we can apply (14.2) 
back to (12.13). As a result we obtain the following equality: 

n n n 



n n n / ^ \ 
r =l q=l ^ 



n n n 



(14.3) 



EEE(^l+^)p, 8 v Tj Mj. 



= 1 r=l 3=1 



This means that the equation (12.13) splits into two parts, first part is (14.3), while 
second part leads to the equations (14.2). 

Both sides of (14.3) do vanish for our special hypersurface S above at its fixed 
point p — po- However, for arbitrary hypersurface S they are nonzero, therefore 
we are to study (14.3) in order to derive other equations for covector field Q. Let's 
denote by B extended tensor field with components 

B * = E E p i (p q % + v q Q>) p °- ( 14 - 4 ) 



k=l 3=1 
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It's clear that tensor field B with components (14.4) is an operator field. Due to 
the presence of P s fc and in (14.4) we have the following equalities: 

B = PoB = BoP. (14.5) 

Relying upon (13.1), (13.2), (13.3) and using (14.4), now we can write (14.3) as 

b(T i ,Br j )=b(T j ,Br i ). (14.6) 

As we noted above, vectors and Tj can be replaced by two arbitrary vectors X 
and Y, see formulas (14.1). Then (14.6) is transformed to 

6(X,B.P(Y))=6(Y,B.P(X)). 

Due to (14.5) and theorem 13.1 we can further simplify this relationship: 

6(X,B(Y))=6(B(X),Y). (14.7) 

Formula (14.7) means that operator B is symmetric with respect to bilinear form 
(13.3). Now we are to utilize this equality. 

Let's denote by W = ImP the image of projection operator P = P(go) for some 
fixed point q n = (po,p) of cotangent bundle T*M (see theorem 13.2). Then W is 
(n — l)-dimensional subspace in T pa {M). It coincides with null-space of covector 
p. Due to (14.5) subspace W is invariant under the action of operator B. More- 
over, operator B is completely determined by its restriction to W. For instance, if 
restriction of B to W is identical operator in W, then B = P: 

B| =id w implies B = P. (14.8) 

Bilinear form b is also completely determined by its restriction to subspace W. 
This follows from relationships (13.4). The equality (14.7) means that restriction 
of operator B to W is symmetric with respect to restriction of b to W. Now 
recall theorem 13.2. It means that arbitrary quadratic form in subspace W can be 
obtained as second fundamental form of some hypersurface S tangent to W. Hence 
restriction of B to W is symmetric with respect to all quadratic forms in W. It 
takes place if and only if the restriction of B to W is a scalar operator, i. e. 

B| = A • idw . 
\w 

This is easy result in linear algebra. Now, applying (14.8) to the above equality, 
for operator B itself we derive the following representation: 

B = AP. (14.9) 

Here A is some scalar factor. It can be expressed explicitly through trace of B: 



(14.10) 
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Formula (14.9) is important result. Combining it with (14.4), we obtain 

E E p i (p q tt + ^ qQk ) p ° = A p; - (14 - n) 

fe=l g=l ^ ' 

Substituting (14.10) for scalar factor A in (14.11), we get the following equality: 



fc=lq=l V " ' 
k=l 0=1 V 



(14.12) 

/n \ pk pr 



This is another additional normality equation. Both (14.2) and (14.12) form a 
system of partial differential equations for components of extended convector field 
Q that defines Newtonian dynamical system in form (7.5) relative to Hamiltonian 
dynamical system with Hamilton function H. Having derived these equations, wc 
proved the following theorem. 

Theorem 14.1. Additional normality condition stated in definition 7.2 for New- 
tonian dynamical system (7.5) in multidimensional case n 3 is equivalent to the 
system of additional normality equations (14.2) and (14.12) that should be fulfilled 
at all points q = (p, p) of cotangent bundle T*M , where p ^ 0. 

15. Connection invariance. 

Deriving additional normality equations (14.2) and (14.12) above we used some 
symmetric extended connection I\ Components of this connection are present in 
(14.2) due to covariant derivatives V r and V s . However, above we did not spec- 
ify which particular symmetric connection T is used. This means that additional 
normality equations (14.2) and (14.12) should be invariant under transformations 

where T* are components of symmetric extended tensor field of type (1, 2). Differ- 
ential equations (14.12) are obviously invariant under transformations (15.1) since 
momentum gradient V is defined without use of connection components (see for- 
mula (10.3)). In general, differential equations (14.2) are not invariant under these 
transformations. However, if component of covector field Q satisfy differential 
equations (14.12), then equations (14.2) become invariant under transformations 
(15.1). In other words, this means that differential equations (14.2) are invariant 
under transformations (15.1) modulo differential equations (14.12). This fact can 
be checked up by direct calculations. 

16. Weak normality equations. 

Let's fix some point qo = (po, Po) of cotangent bundle T*M with p ^ 0. It yields 
initial data for Newtonian dynamical system written in form (7.5) and defines a 
trajectory p = p(t) of this dynamical system passing through the point p = pa 
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at time instant t — 0. Null-space of covector p is a hyperplane in tangent space 
T po (M). Let's denote it by a. One can draw various hypersurfaces passing through 
the point po and tangent to hyperplane a at this point. Suppose that S is one of 
such hypersurfaces and suppose that n = n(p) is smooth normal covector of S in 
some neighborhood of the point po- At the very point p = po we have the equality 

Po = vq ■ n(p ), where v Q ^ 0. (16-1) 

Now let's take some smooth function on S normalized by the condition (7.15) and 
set up initial data (7.4) for Newtonian dynamical system (7.5). Solving Cauchy 
problem with these initial data, we get a family of trajectories for dynamical system 
(7.5), which includes our initial trajectory passing through the point p = p a . This 
is easily seen if we compare (16.1) and (7.15). 

In local coordinates the family of trajectories constructed just above is repre- 
sented by functions (2.4). Using them, we define variation vectors (2.5) and devia- 
tion functions (2.6). Now suppose that Newtonian dynamical system (7.5) satisfy 
strong normality condition (see definition 8.1). This means that by proper choice 
of function v{p) we can make all deviation functions ipi, ... , f n -i to be identically 
zero. Hence initial conditions (7.2) are fulfilled. From section 7 we know that ini- 
tial data (7.2) are equivalent to Pfaff equations (7.13) for v. We can vary constant 
iso ^ in normalizing condition (7.15), and for each value of this constant due to 
strong normality condition we would have some function v{p) on S satisfying Pfaff 
equations (7.13). Due to lemma 7.2 then Pfaff equations (7.13) are compatible in 
the sense of definition 7.1. Thus we have proved the following theorem. 

Theorem 16.1. Strong normality condition implies additional normality condition 
for Newtonian dynamical system (7.5). 

Additional normality condition is formulated in definition 7.2. In section 14 we 
have shown that additional normality condition is equivalent to additional nor- 
mality equations (14.2) and (14.12). Thus strong normality condition leads to the 
equations (14.2) and (14.12) for components of covector field Q. However, it can 
yield much more. Indeed, it implies initial condition 

(f>i\ =0 (16.2) 
lt=o 

in addition to initial conditions (7.5). Let's calculate second derivatives ipi for 
deviation functions <pi, ... , tp n -i by differentiating formula (7.10). First of all 
note that formula (7.10) itself can be written in terms of covariant derivatives 

* = - E ¥ -E(^r- q) <■ (i6.3) 

s—1 s—1 ^ ' 

Note that the equations of Newtonian dynamics in form (7.5) also admit covariant 
derivatives instead of partial derivatives in them: 

i s = ^, VtPs = -^- + Q s . (16.4) 
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Applying covariant derivative V Xi to (16.4) and denoting V^p = £j, we obtain the 
following differential equations for components of vector n and covector : 

^? - E ^i^r) ^ + E v {^) ( 16 - 5 ) 



lb It IL s ^ 

r— 1 o—l m—1 ^ 



n n 



+ E E E Kr < + E = (16-6) 

g— 1 m— 1 \ r— 1 r— 1 / 

r=l v 7 r=l v 7 

Here D™ 9 and D™ r are components of dynamic curvature tensor D given by formula 

dT k 

D% = -^L. (16.7) 

Tensor D has no analogs in Ricmannian geometry since its components (16.7) do 
vanish for non-extended connections. In (16.6) we have quantities R™ qr which are 
components of another curvature tensor R given by formula 

rij ~ Q x i Q x j 2-^1 ™ jr 2-^1 3 m ir 

m—1 m—1 

(16.8) 

n n p, v k n n „ pfe 

+ W v r a — - W v r a -i^ 

m—1 a— 1 m—1 a— 1 

In Riemannian geometry (16.8) reduces to standard formula for components of 
Riemann curvature tensor. 

Now let's differentiate (16.3) with respect to time variable t. It is equivalent to 
applying covariant derivative V t to to this equality: 

- 1 V ^ - 1 (16.9) 

-Ev t (^-Qf/-E(^-Q s )v./. 

s— 1 x 7 s— 1 v 7 

Substituting (16.5) and (16.6) into (16.9), we obtain the following equality for iff. 



n / n 



r—1 \ s—1 s—1 
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" V S H ~„ V r O A / A V,Q V S H 



r=l s=l \ / r=l \ s=l 

71 X7 S Q f V H \\ n n , 

^~^{ — h~ + Qs ) v ^ <+££( v ^-v r Q s 

s— 1 ^ ' / r— 1 s— 1 ^ 



r— 1 s— 1 x 



Terms with curvature tensors are canceled due to the following identities: 



[Vi, Vj]H = - £ £?>fc k % V S H, (16.10) 

fe=l S = l 

n n 

[v if v^iH = £ i> ^ ™ (i6.il) 



fe=l s=l 



The identities similar to (16.10) and (16.11) in v-represcntation were derived in 
Chapter III of thesis [6]. 

Thus, formula for is derived (see above). It is rather huge. In order to simplify 
this formula we introduce two extended fields a and (3 with components 



a 



n 2 n n 2 V n 

s=l s=l v 



E— ( v ^ +— 



(16.12) 



/ A v a n v s # A v s ^ / v s i? ^ » , „ „ 



ft 2 n ^ V n 

= 1 s=l v 



(16.13) 



+ £(V s Q r -V r Q s - ^Qs) ^--E(^T-^j V 

s— 1 ^ / s _i ^ ' 

Using notations (16.12) and (16.13), we can write formula for (pi as follows: 

n n 

^ = ]>> s V n p s + £/3 s T/. (16.14) 

s=l s=l 

Then, using projector P with components (12.3), we can transform this expression: 

n n n n 

& = £E arp 'v T( p,+E;E>i?T/+ 

r— 1 s— 1 r— 1 s=l 

(16.15) 

+E E aT vr v Ti p, +EE^ir p « ^ 

r— 1 s— 1 r— 1 s — 1 
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Comparing (16.15) with formula (16.3) for ipi, we can write formula (16.15) as 

n n n n 

r—1 s—1 r—1 s—1 

(16.16) 

-i>i«> (^-^)<+ex>^w 

s=l ^ ' r=l s=l 

Now let's introduce other two extended fields a and rj with components 

»7r=/?r-<P|a>(^-Qr), CT = E ^ ( lfU7 ) 

Then, taking into account formula (2.6) for t^, we can write (16.16) as follows: 

n n n n 

^ + (p\ a )&-a V , i = j2Yl aT p r v ^p s + EE ^ p s T i- ( 16 - 18 ) 

r—1 s—1 r—1 s—1 

Combining (16.18) with (7.2) and (16.2), we obtain the equality 

n n n n 

E E " r p r v ^ Ps + E E * p s T ° = °- ( 16 - 19 ) 

r—1 s — 1 r—1 s — 1 

If we remember operator b determined by (13.1) and (13.2) and if we use formulas 
(13.3) and (13.4) for bilinear form b, then (16.19) is written as 

- (b(Ti) | Pa) + (t) I Pt,) = 0. (16.20) 

Now recall that in the beginning of this section we have taken a point go = (po, Po) 
and considered a set of hypersurfaces in M passing through the point po tangent 
to null-space of covector po ^ 0. Therefore, when equality (16.20) is written for 
the point po, we can replace Pii by PX, where X is an arbitrary vector of tangent 
space T po (M). In a similar way, due to theorem 13.2, covector b(Tj) in (16.20) can 
be replaced by arbitrary covector y e T* g (M). Hence (16.20) breaks into two parts 

(y|Pa)=0, <77|PX>=0 (16.21) 

with arbitrary vector X and arbitrary covector y. In coordinate form these two 
equalities (16.21) are equivalent to the following ones: 

n n 

E P > S=0 > J2 P r r >s=0- (16.22) 

8=1 8=1 



Looking at (16.12), (16.13), and (16.17), we see that the equalities (16.22) form 
a system of partial differential equations for components of covector Q written in 
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terms of covariant derivatives V and V. They are called weak normality equations. 
Let's write them explicitly. For the first equation (16.22) we have 



n n 

EE 



Y7STT f \7 r O \ 

* " 'v r Q s + —Q s )P r q = 0. (16.23) 



r—l 5=1 \ 



Second equality (16.22) leads to more huge equations 

EE V s Q r + ^-Q r -V r g s + ^-Q s + 

r=l s= l \ V " - / 




(16.24) 



Theorem 16.2. Strong normality condition for Newtonian dynamical system 

(16.3) implies weak normality equations (16.23) and (16.24) to be fulfilled at all 
points q = (p,p) of cotangent bundle T*M , where p / 0. 

Note that theorem 16.2 is stated for Newtonian dynamical system written in form 

(16.4) . However, weak normality equations (16.23) and (16.24) are invariant under 
the transformations (15.1) (this can be proved by direct calculations). Therefore 
theorem 16.2 is equally applicable to Newtonian dynamical system written in form 
of the equations (7.5). 

17. Equivalence of strong and complete normality conditions. 

In section 16 we have derived weak normality equations (16.23) and (16.24) from 
strong normality condition for Newtonian dynamical system (see definition 8.1). 
Here we reveal their relation to weak normality condition considered in section 6. 
As in section 6, let's consider one-parametric family of trajectories p = p(t, y) of 
Newtonian dynamical system (16.4). In local coordinates it is represented by func- 
tions (6.1). Differentiating them with respect to parameter y, we obtain variation 
vector r, see formula (6.3). Then we can define deviation function (6.4). Unlike 
section 6, in the above calculations in section 16 we used p-representation rather 
than v-representation. Indeed, the equations of Newtonian dynamics (16.4) and 
all normality equations (14.2), (14.12), (16.23), and (16.24) are written in terms 
of Hamilton function H and in terms of covariant derivatives (10.3) and (10.7), 
which require momentum representation of extended tensor fields. For this reason, 
instead of functions l (t,y) in (6.5), we consider components of vector £ = \7 T p: 

c\ n n 



k=l 9=1 



(compare with formula (12.9) above). Functions r 1 , ... , r", £ x , . . . , £™ considered 
as functions of time variable t for fixed y satisfy a system ordinary differential 
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equations which are quite the same as the equations (16.5) and (16.6) above, except 
for the absence of index i now: 

V-- = £v^k+£vY™U (17.1) 



r=l 



-Q q )r r + 



n n n / T7 if 

r—1 q—1 m—1 ^ 

+EE E Kr + E D 7s r = (17.2) 

q—1 m—1 \ r—1 r—1 / 

-t*<¥-«-)fr-t^¥-«-)''- 

These equations (17.1) and (17.2) can be understood as linearizations of the equa- 
tions of Newtonian dynamics (16.4). They are linear equations with respect to 
functions r 1 , . . . , r™, . . . , but with non-constant coefficients. Their coef- 
ficients are functions of time variable t determined by a trajectory p = p(t) of 
Newtonian dynamical system (16.4). For a fixed trajectory p = p(t), i.e. when 
parameter y in p = p(t, y) is fixed, solutions of the equations (17.1) and (17.2) form 
n-dimensional linear space. We denote this space by 1. In essential, it is the same 
space T as in section 6. 

Let's consider deviation function ip determined by formula (6.4) and its time 
derivatives. Now ip and functions r 1 , . . . , r™, . . . , £™ are not related to any 
hypersurface S. Nevertheless, repeating the same steps as in deriving formula 
(16.18), one can derive the following equality for deviation function ip: 

n n n n 

y5+(p\ a )^-a^ = J2J2 arp r^+J2J2^Plr s . (17.3) 

r—1 s—1 r—1 s—1 

Remember that weak normality equations (16.23) and (16.24) are expanded form 
of the equations (16.22). Therefore, if weak normality equations are fulfilled, then 
(17.3) reduces to second order ordinary differential equation for ip: 

Cp+ (p| a) ip - aip = 0. (17.4) 

Comparing (17.4) with (6.11), we find that Newtonian dynamical system fits the 
definition 6.1 under the condition that weak normality equations (16.23) and (16.24) 
for extended covector field Q arc fulfilled. 

Converse result is also valid, i. e. weak normality condition stated in definition 6.1 
implies weak normality equations (16.23) and (16.24) to be fulfilled. Let's prove 
it. Suppose that we take some trajectory of p = p(t) of Newtonian dynamical 
system (16.4). This determines coefficients in linear differential equations (17.1) 
and (17.2) for components of r and £ and fixes linear space T. Then, according 
to definition 6.1, for each solution of the equations (17.1) and (17.2) corresponding 
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deviation function if should satisfy differential equation (6.11). For time derivatives 
if and ip we have the following equalities: 

S— 1 5=1 X 7 



E« s ?+E^ tS - ( 17 - 6 ) 



(compare with (16.3) and (16.14) above). For the function f itself we have formula 
(6.4). Due to differential equation (6.11) from (17.5) and (17.6) we derive 

According to definition 6.1, this equality should be fulfilled for all solutions of 
differential equations (17.1) and (17.2). Therefore 

V S H 

a s +A—=0, (17.7) 
(3s + a(^--q)\ -B Ps =0. (17.8) 



Multiplying (17.7) by p s and summing over s running from 1 to n, we obtain 



A=-J2a s Ps = -(p\a). (17.9) 

s=l 

In a similar way, multiplying (17.8) by V S H and summing over s running from 1 to 
n, we derive formula for coefficient B in (6.11): 

»- ± (* + ^ (¥-«•))-'■ (17 - 10 » 

Now, let's substitute (17.10) and (17.9) back into the equations (17.7) and (17.8). 
Then let's multiply (17.7) and (17.8) by P s 9 and respectively and sum them over 
index s. This yields the equalities (16.22) which are equivalent to weak normality 
equations (16.23) and (16.24). Thus we have proved the following theorem. 

Theorem 17.1. Weak normality condition stated in definition 6.1, when applied 
to Newtonian dynamical system written as (7.5), is equivalent to the system of 
weak normality equations (16.23) and (16.24) that should be fulfilled at all points 
<1 = (p,p) °f cotangent bundle T*M , where p^O. 

Combining theorems 16.1, 16.2 and 17.1, we obtain another theorem. 

Theorem 17.2. Strong and complete normality conditions are equivalent to each 
other either for n = 2 and in multidimensional case for n ^ 3. 
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18. Summary and conclusions. 

Primary goal of present paper is to generalize theory of Newtonian dynamical 
systems admitting normal shift from Riemannian geometry to the geometry de- 
termined by some Lagrangian or, equivalently, by some Hamiltonian dynamical 
system. In the above sections this goal is reached in essential: 

- we have found proper statement for the concept of normal shift in Lagrangian 
geometry and have defined class of Newtonian dynamical systems admitting 
normal shift of hypersurfaces; 

- we have derived complete system of normality equations thus obtaining effective 
tool for studying this class of dynamical systems. 

However, some details of constructed theory appear to be different from those 
one could expect in the beginning. Indeed, being generalization of Riemannian 
geometry, geometry of Lagrangian dynamical system could have some connection 
canonically associated with Lagrange function. But even if such connection does 
exist, theory of normal shift does not reveal it. All normality equations (14.2), 
(14.12), (16.23), and (16.24) are invariant under transformations (15.1). Therefore 
they are of connection-free nature. There is a problem of writing them in coordinate 
covariant tensorial form without use of connection at all. This problem will be 
considered in separate paper. 
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